Observation of the Decay K~ 



tv n + n and Measurements of the Branching Ratios 

for K — > 7T ± /_i + /l _ 
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Using data collected with the HyperCP (E871) spectrometer during the 1997 fixed-target run 
at Fermilab, we report the first observation of the decay K~ —* ir~ n + fi~ and new measurements 
of the branching ratios for K — > tt /j, + fj,~. By combining the branching ratios for the decays 



n+fi+Li- and K~ -> n~ y+ , we measure Y{K ± -> 7r ± ^ + /i-)/r(A" ± -> all) = (9.8 ± 
8 . The CP asymmetry between the rates of the two decay modes is [r(_K" + — 1 



K+ 

1.0 ± 0.5)xl0 
■K+H+yT) -T(K 



-0.02 ±0.11 ±0.04. 



PACS numbers: 13.20.Eb, 14.40.Aq, 11.30.Er 



The rare decay of charged K mesons to a pion and a 
lepton pair (K vll , where I = e or fx) can be used to study 
flavor-changing neutral currents as a higher-order process 
in the standard model and to explore new physics. How- 
ever, to realize these goals, it is necessary to understand 
the dominant K — > 71^7* radiative transition which 
involves long-distance hadronic effects. Since it is diffi- 
cult to calculate such effects, chiral perturbation theory 
(ChPT) including electroweak interactions has been ap- 



plied to K^ u [|1[ . In a recent model-independent analysis 
with ChPT at 0{p 6 ) §J, the K+ n decay rate and the 
form factor characterizing the dilepton invariant-mass 
spectrum are calculated in terms of two parameters a+ 
and . Determining these parameters from the mea- 
sured branching ratio and the dilepton mass spectrum 
of the ifjJee decay, this analysis predicts the ratio 1Z = 
B(K+ flfi )/B(K+ ee ) to be greater than 0.23 and an in- 
crease in the CP asymmetry between the decay rates of 
K^ u and K~ u compared to the leading-order estimation 
in the chiral expansion . The resulting expected CP 
asymmetry is ~ 10~ 5 [§,]3|. 

The Kiee decay has been studied by several exper- 
iments H |, §, and the Particle Data Group (PDG) 
has compiled a mean branching ratio B(K+ ee ) — (2.88 ± 
0.13) x 10~ 7 {§. The decay was first observed 

by the E787 Collaboration at the Brookhaven National 
Laboratory (BNL). With 13 fully reconstructed 3-track 



and 196 partially reconstructed 2-track events, they de- 
termined B{K+ ) = [5.0 ± 0.4(stat) ± 0.7(syst) ± 
0.6(theor)] xlO" 8 @. The E865 Collaboration at BNL 
has subsequently observed 430 fully reconstructed K+„„ 
events and measured B(K+^) = [9.22 ± 0.60(stat) ± 
0.49(syst)] xl0~ 8 The discrepancy between these 
two experimental results is more than three standard de- 
viations and is not understood. In addition the value of 
1Z obtained with the E787 result and the average value of 
B(K+ ee ) is inconsistent with the prediction and is diffi- 
cult to accommodate within the standard model j^j . It is 
thus important to resolve the discrepancy in the B(K i, „) 
measurements . 

In this Letter we present results on the first observa- 
tion of K~ — » 7r~/i + ^~ decay, the branching ratios of 
^nuu decays, and a limit on their CP asymmetry. The 
HyperCP (E871) experiment in the Meson Center beam 
line of Fermilab, which was designed primarily to study 
CP violation in S~/S + and A/A decays, took data in 
1997 and 1999. We report here on an analysis of the 
1997 data. A plan view of the spectrometer is shown in 
Fig. 0; a detailed description is given elsewhere JTo| . 

K ±: s were produced by directing an 800-GeV/c pro- 
ton beam, with a typical intensity of 1.2 x 10 11 protons 
per 20-second spill, onto a 2-mm-square copper target. A 
target length of 22 (60) mm was used when taking pos- 
itive (negative)-secondary-beam data. The secondary- 
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FIG. 1: Plan view of the HyperCP spectrometer. The length 
(width) of the spectrometer is about 62 m (4 m). 



beam momentum was selected by a curved collimator 
embedded in a 6-m-long dipole magnet (Hyperon mag- 
net) with a horizontal magnetic field of 1.67 T. The mean 
momentum of the secondaries was about 170 GeV/c, and 
the charge of the secondary beam was selected with the 
sign of the Hyperon magnet field. The typical secondary- 
beam rate was 13 MHz at the exit of the collimator. Fol- 
lowing a 13-m evacuated decay pipe were eight multi- 
wire proportional chambers (C1-C8) with wire spacings 
increasing from 1 to 2 mm. A dipole magnet (Analyz- 
ing magnet) deflected charged particles horizontally with 
a transverse-momentum kick of 1.43 GeV/c. Particles 
with the same charge as the secondary beam were de- 
flected to the left and those with opposite charge to the 
right. At the rear of the spectrometer were two scintil- 
lation hodoscopes and a calorimeter, used for triggering, 
and a muon detector system. The muon system consisted 
of two similar detectors on either side of the secondary 
beam line. The muon detectors had three layers of 0.81- 
m-thick steel absorber, each layer followed by a track- 
ing station consisting of vertical and horizontal planes 
of proportional tubes with a pitch of 25.4 mm. At the 
rear of the third layer, behind the proportional tubes, 
were two orthogonal scintillator hodoscopes used for the 
muon triggers. 

The rates were measured relative to the K ± — > 

7r ± 7r + 7r _ (K^ 3 ) rates (the normalization mode). The 
trigger for the K^ 3 events was a heavily-prescaled coinci- 
dence of at least one hit in each of the left- and right-side 
hodoscopes. The trigger for the events, the unlikc- 

sign dimuon trigger, required the presence of at least one 
hit in each hodoscope in the muon detectors in coinci- 
dence with hits in the left- and right-side hodoscopes. 
The sign of the secondary beam was changed every few 
hours by reversing the polarities of both the Hyperon and 
Analyzing magnets. The typical run cycle was two runs 
with positive secondary beam followed by one run with 
negative secondary beam. 



Data for the and K n3 triggers were processed 

with the same reconstruction program. A three-track 
event topology was required. The sample of ifi^ can- 
didates was selected by requiring a muon track in each 
muon station and an unlike-sign dimuon trigger. The 
muon track was required to have hits in at least two of 
three muon chambers in both x and y views with corre- 
sponding in-time hits in the muon hodoscope. The al- 
gorithm for muon selection was checked using muon cal- 
ibration data. The momentum of accepted muons was 
greater than 20 GeV/c. 

Similar selection criteria were applied to both the sig- 
nal (-K^uu) and normalization (K^ 3 ) event samples. One 
right-side track and two left-side tracks were required. 
Two pions from the K^ 3 decays were required to project 
to the same fiducial areas in the left- and right-side muon 
hodoscopes as the muons from the decays. The to- 

tal momentum of the reconstructed tracks was required 
to lie between 120 and 250 GeV/c, consistent with the 
momentum spread of the secondary beam. To ensure 
that the kaons came from the target, the combined three- 
track momentum had to point back to within 5 mm of the 
center of the target in both x and y. The decay vertex of 
the three tracks was calculated by the method of distance 
of closest approach using the front track segments, and 
was required to be well within the vacuum decay region. 

The and K^ 3 event candidates were required to 

have a decay topology consistent with a single vertex. 
This was done by first fitting the three tracks to a sin- 
gle vertex using the wire hits in C1-C4 and computing 
the x 2 of the fit. In addition the average separation be- 
tween the three pairs of tracks in the x — y plane at the 
z position of the vertex determined by distance of closest 
approach was measured. Clean K^ 3 decays were used to 
study selection cuts based on the vertex-constrained-fit 
X 2 and the average separation distance. Candidates were 
accepted if x^/ndf < 2.5 and the average separation dis- 
tance was less than 2 mm. These vertex cuts provided 
strong rejection of backgrounds, in particular those from 
two- vertex hyperon decays. 

For events passing these cuts, the reconstructed three- 
track invariant mass M vflfl distributions with a de- 
cay hypothesis are shown in Fig. ||. The signals in the 
vicinity of the kaon mass are clear and unambiguous. 



To determine 

r(ir± 

was used: 
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the branching ratios B(K* ) = 
") /r(ii' ± — > all) the following formula 
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where the N° bs, s are the numbers of observed events, 
the Ai's are the geometrical acceptances, the ej's are 
the event-selection efficiencies, e'" ; s is the relative trigger 
efficiency, and 200 is the prescale factor. The branch- 
ing ratios B(K± 3 ) are (5.59 ± 0.05)% @. Unbinncd 
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FIG. 2: The 7r/x/i invariant mass distributions for (a) positive- 
and (b) negative-secondary-beam data after cuts. The large 
peak comes from K^ 3 events with in-fiight decays of pions and 
pion punch-through to the muon stations. The inset shows 
the region around the kaon mass in more detail; points are 
data and the solid lines show the results of the likelihood fits. 



maximum-likelihood 

470 MeV/c 2 < 



fits for the events in the region 
< 520 MeV/c 2 were used to esti- 



mate the numbers of events, where the mass range 

for the fit was chosen to minimize the statistical error. 
Gaussian and linear functions were used to fit the sig- 
nal and background respectively. The numbers of sig- 
nal events from the fits were iV° b + = 65.3 ± 8.2 and 

jyobs 



= 35.2 ± 6.6, where the errors are statistical. Af- 



ter background subtraction, the observed numbers of nor- 
malization events were A° fc f = (4.446 ± 0.010) xlO 5 and 

K ir3 

N°J?1 = (2.318 ±0.008) xlO 5 , where the errors are domi- 
nated by the background estimation. 

Background contributions to the samples from 

pion punch-through in the muon stations and in-flight de- 
cays of pions in K^ 3 decays and from — > tt + it~ p^v^ 
were negligible, mainly due to the good resolution of the 
spectrometer. Monte Carlo (MC) studies for these de- 



cays showed that the reconstructed masses with the i:\i\x 
hypothesis of these modes were well below 480 MeV/c 2 . 

Many checks to verify the MC simulation were made. 
The kinematic parameters for K ± production at the tar- 
get were tuned separately for the K^ 3 and K~ 3 normal- 
ization samples. The PDG Dalitz parameters for K^ 3 
decays were used in the MC program to properly sim- 
ulate this decay. The distributions for kaon momenta 
and vertex positions of K^ 3 decays inside the decay pipe, 
hit positions of pions from the decays, and reconstructed 
masses were compared with data. They all showed good 
agreement. 

In the simulation of decays, we assumed a vec- 

tor interaction, based on the results of previous exper- 
iments [[| [|, and so the angular distribution of ^+ in 
the rest frame of the virtual photon was generated ac- 
cording to the expression (m 2 ±p 2 sin 2 9), where m M and 
are the mass and momentum of the muon, and 9 is 
the angle between the directions of the fi + and the pion 
in this frame. In addition, we employed a form factor 
as described below. The dilepton invariant mass (My) 
spectrum from the K^ u decays is described by the form 
factor </>(Mj 2 ) = (j) (l + 8MfJm 2 K+ ) where S and <j) are a 
slope parameter and a constant, respectively, and m^t 
is the mass of the charged kaon || [fjj. To extract <5, 
we averaged the measured values of the slope parameters 
determined by BNL E777 using K+ ee decays || and by 
BNL E865 with K+ ee [| and events §. We ob- 

tained S — 2.08 ±0.12, where the error is only statistical. 
Figure |^ compares data for the dimuon invariant mass 
for decays with simulation results for 3 values of 5. 

From the MC simulation the geometric acceptance for 
K^ 3 and K^- decays inside the vacuum decay region for 
positive (negative)-secondary-beam data were estimated 
to be 94.4% (94.2%) and 47.5% (47.7%), respectively, 
where we have taken into account pion decay in flight 
and pion punch-through to the muon stations. Analyses 
identical to those applied to the data were used for the 
simulated normalization and signal events that passed 
the trigger requirements. The event selection efficiencies 
for K± 3 and decays were 77.9% (76.0%) and 80.3% 

(78.3%) for positive (negative) secondary beam data re- 
spectively. The difference in efficiencies for K + and K~ 
events is primarily due to the slightly different momen- 
tum spectra for K + and K~ production. To estimate 
the trigger and dimuon selection efficiencies, the full data 
sample was used. Since the trigger for the normalization 
mode was common to that for the signal mode, we could 
determine the relative trigger efficiency for the signal 
mode with respect to the normalization mode. The rela- 
tive trigger efficiency e'" ; s was found to be (86.6 ± 2.7)%. 
The dimuon selection efficiency e^+^- was (93.8 ±0.3)%. 

A summary of the systematic uncertainties is given 
in Table Systematic uncertainties on the measured 
branching ratios were minimized by using the same cuts 
both for the signal and for the normalization modes, with 
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FIG. 3: The invariant mass after background subtrac- 

tion for (a) positive- and (b) negative-secondary-beam data 
after cuts. Points are data. The solid, dashed and dotted 
lines represent the simulation result with S — 2.08, 1.72 and 
2.44, respectively. 



the exception of the muon-identification cuts which were 
applied to the signal events only. The major systematic 
uncertainties were the background estimation for the sig- 
nal events and the trigger efficiency. In order to estimate 
the uncertainties in the numbers of signal events, various 
mass ranges were used for the fits, and fits with con- 
stant and quadratic functions for the background were 
used. The variations of the beam position at the tar- 
get and the magnetic field strength in the Hyperon and 
Analyzing magnets during the run were monitored. By 
varying the assumed beam position, magnetic fields, K 
production parameters, and the decay parameters in the 
MC, we determined the sensitivity to these parameters of 
the relative acceptances and efficiencies for the signal and 
normalizing modes. Realistic ranges for these parameters 
were then used in computing the MC-dependent system- 
atic uncertainties in Table |. Finally, the effect of the 
uncertainty in the slope parameter 8 was estimated by 
varying it by three standard deviations (see Fig. ^), and 
found to be negligible. 

Based on Eq. (|l|) we determined the branching ratios 
to be B(K+^) = (9.7±I.2±0.4) x I(T 8 and B{R-^) = 
(10.0±L9±0.7) x I0~ 8 , where the first and second errors 
are statistical and systematic respectively. In terms of 
the partial decay rates ), the CP asymmetry in 

these modes is defined as 



T{K+) - T(K-) 



(2) 



From our measurements of B(K^ ), we obtained 



A(K^) = -0.02 ± O.ll(stat) ± 0.04(syst). 

In conclusion, we have observed both K + — > 7r + /i + /i~ 
and K~ — > 7r decays and measured their 

branching ratios. This is the first observation of the 
K~ — > 7T - /i + [i~ decay. The CP asymmetry parameter, 
A(2f^„„), has been extracted from these measurements. 
Our result is consistent with no CP violation. Assum- 
ing CP symmetry is valid, we can combine B(K+ ) 
and B(K~^^) to give a result of [9.8 ± l.O(stat) ± 
0.5(syst)] x 10 -8 , which is 3.2 standard deviations higher 
than the BNL E787 measurement || , but consistent with 
the BNL E865 result § and, using the PDG value for 
B(K+ ee ), is also consistent with the model-independent 
analysis for the 1Z value S . 



Source 


a B /B (%) o B 


(%) 


Beam targeting 


1.1 


(0.9) 


1.0 


Magnetic field 


0.5 


(0.4) 


0.5 


Trigger efficiency 


3.1 


(3.1) 


3.1 


Muon identification 


0.3 


(0.3) 


0.3 


Background Estimation (K^r 3 ) 


0.2 


(0.3) 


0.2 


Background Estimation (K^pp) 


2.4 


(6.4) 


3.6 


Data and MC disagreement 


0.4 


(0.4) 


0.4 


Dalitz parameter 


0.3 


(0.3) 


0.3 


Slope parameter (S) 


0.2 


(0.3) 


0.2 


B(K+ -> tt+tt+tO 


0.9 


(0.9) 


0.9 


Total 


4.2 (7.3) 


5.0 



TABLE I: Summary of estimated systematic uncertainties 
in the branching ratios of B(K^^), and, in parentheses, 
B{K^^^). The last column gives the systematic errors for 
the combined branching ratio, B(K^„ I1 ). 
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